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Abstract 

We exploit mappings between quantum and classical systems in or¬ 
der to obtain a class of two-dimensional classical systems with critical 
properties equivalent to those of the class of one-dimensional quan¬ 
tum systems discussed in a companion paper [6]. In particular, we use 
three approaches: the Trotter-Suzuki mapping; the method of coherent 
states; and a calculation based on commuting the quantum Hamilto¬ 
nian with the transfer matrix of a classical system. This enables us to 
establish universality of certain critical phenomena by extension from 
the results in [6] for the classical systems identified. 


1 Introduction 

Mappings between statistical mechanical models have provided new pathways 
to compute thermodynamic properties of systems which were previously in¬ 
tractable mmm- In particular, critical phenomena in d-dimensional quan¬ 
tum systems have been investigated by mapping them to (d+ l)-dimensional 
classical systems for which there are better developed techniques, such as 
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Monte Carlo simulations [10J 16]. For example, one well known connection is 
that between the one-dimensional XYZ model and the two-dimensional zero- 
held eight-vertex model, namely that the Hamiltonian of the quantum model 
and the transfer matrix of the classical model have the same eigenvectors. 
Baxter [1] found the ground state energy for the XYZ model by first finding 
the partition function of the eight-vertex model and then showing that the 
quantum Hamiltonian is effectively the logarithmic derivative of the transfer 
matrix for the classical system. 

In this paper we exploit these quantum to classical (QC) mappings for the 
opposite reason: to take advantage of known ground state critical behaviour 
in a general class of quantum spin chains to determine the finite-temperature 
critical properties of an equivalent class of classical spin systems. 

In [ 6 ] we computed the critical exponents s, v and z, corresponding to the 
energy gap, correlation length and dynamic exponent respectively, for a class 
of quantum spin chains, establishing universality for this class of systems. 
We also computed the ground state correlators (crfcrf +r ) , (ofcrf +r ) and 
(n,:=i a i) g f° r this class of systems when translation invariance is imposed. 
These correlators were found to exhibit quasi-long range order behaviour 
when the systems are gapless, with a critical exponent dependent upon the 
system parameters. 

The class of quantum spin chains studied in [ 6 ] consists of M spin-1/2 
particles in an external field h, with a Hamiltonian quadratic in Fermi oper¬ 
ators given by 

M M 

H = Y,( A s$ bt+ l B #(b%-l>jb t ))-2hY i b]b j , ( 1 . 1 ) 

j,k j =1 

where the bj s are the Fermi operators satisfying the usual Fermi commutation 
relations 

{&;,!>*} = in,, {j{ 4 } = {b„h} = (f >») 2 = (b,f = 0 . ( 1 . 2 ) 

The measure of anisotropy 7 is real, with 0 < 7 < 1; the matrix A ]k must be 
Hermitian and Bj k antisymmetric, both containing only real entries without 
loss of generality; and periodic boundary conditions &M+j = bj are assumed. 

This model can be diagonalised HH so that 

^ = ^2\ A q\ rfqVq + C, (!-3) 

q 

with the dispersion relation |A 9 | determined by matrices A ]k and B 3kl the 
r] q s are Fermi operators, and C is a constant. 
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Iii 00 Keating and Mezzadri restricted the Hamiltonian (1.1) to pos¬ 
sess symmetries corresponding to the Haar measure of each of the classical 
compact groups U(N), 0 + (2N ), Sp(2N), 0 + (2N + 1), 0~{2N + 1) and 
0~(2N + 2), enabling the calculation of \A q \ using techniques from random 
matrix theory. This corresponds to a symmetry classification of spin chains 
similar to that introduced for disordered systems by Atland and Zirnbauer 
cm [19] . These symmetry properties were encoded into the structure of 
the matrices A jk and B jkl as summarised in Table [l] in Appendix [A] For 
example, when restricted to U(N) symmetry^ 0 [9] 


Aq = 4 I T + ^ (a(k) cos kq + ib(k) sin kq ) 

\ k=1 

= 4 (a q + ib q ), 

with real and imaginary parts given by 


(1.4) 


Ctq — 

L 

T + a{k ) cos kq , 

k= 1 

L 

and b q = b{k) sin kq , 

k =1 

(1.5) 

and 

r=* 

2 

fa(0), L = ^, 
\a(0) + (—l)'a(f), 

if M is odd, 

L = f - 1, if M is even, 

(1.6) 

where q is the wave number 




27 tI 

q = —— with 

H M 

l = 0 ,...,M- 1. 

(1.7) 


In general, the symmetry constraints were achieved using real functions 
a(j) and b(j), even and odd functions of Z/MZ respectively, to dictate the 
entries of matrices A ] k and Bj kl as reported in Table [T] in Appendix [A] 

Exploiting the formalism developed in 0 [9] enabled us to compute the 
critical properties of this class of spin chains [6], demonstrating a dependence 
of the critical exponents on system symmetries and establishing universality 
for this class of quantum systems. 

Having established universality for the above class of quantum spin chains 
in 0, we now make use of QC mappings to obtain a class of classical systems 
with equivalent critical properties, establishing universality for this class of 
classical systems as well by extension. This is our main goal. 

A quantum and a classical system are equivalent if their partition func¬ 
tions are the same; such a correspondence, however, is not unique as different 

1 For the other symmetry classes see [8|. 


3 





classical systems can be equivalent to the same quantum system. We shall 
here adopt the following different approaches to map the partition func¬ 
tions of the quantum spin chains (1.1) onto those of a general class of two- 
dimensional classical systems: 


• the Trotter-Suzuki formula (Section [2j; 

• the method of coherent states (Section [3]); 

• the simultaneous diagonalisation of the quantum Hamiltonian and the 
transfer matrix for the classical system (Section [4]). 


2 Trotter-Suzuki mapping 


This approach was developed by Suzuki HZ], who applied the Trotter product 
formula 


e A+B = i e"e" 


A B 


o, 


( 2 , 1 ) 


to map the partition function for a d -dimensional quantum system to that for 
a (d+ l)-dimensional classical one. In particular he applied it to the partition 
function of a d -dimensional quantum Ising model in a transverse magnetic 
held, mapping it to that of a (d + l)-dimensional classical Ising model [l?j. 
He then proved the equivalence of the critical properties of the ground state 
of the quantum system and the finite temperature properties of the classical 
system. 

Here we harness this technique to supply us with a class of two-dimensional 
classical systems with critical properties equivalent to those of the ground 
state of the quantum spin chains (1.1). Like the original quantum system, 


the classical counterparts are also able to possess symmetries reflected by 
those of the Haar measure of each of the different classical compact group^J 
enabling the dependence of critical properties on system symmetries to be 
observed. 

There are many ways to apply the Trotter-Suzuki mapping to the parti¬ 


tion function for the class of quantum spin chains (1.1), resulting in different 


2 This is observed through the structure of matrices A,fc and summarised in Table [I] 

inherited by the classical systems. 


4 








classical partition functions. Those that we obtain are of the form 


all states 

Z £, = 'y ^ g — 0cl'Hcl({ s i,j}) 

restricted states 

z c = yi n»< 

all configurations % 

Pc\Hc\ }) 


Z D - Y1 


(2.2a) 
(2.2b) 
(2.2c) 
(2.2d) 


where T-[ cl is the effective classical Hamiltonian. In (2.2a) and (2.2b) TL 


a real function of the classical spin variables Sij = ±1 and in (2.2d) it is 


is 


a complex function of the classical spin variables = ±1, which 

represent the eigenvalues of the Pauli matrices erf, of and of respectively. 
The function / ({s«,j}) is also a real function of the classical spin variables 


Sij = ±1, and we find that if /({sjj}) = 1, then (2.2a) has the familiar 


form of a classical partition function, with TL C i representing the Hamiltonian 
describing the effective classical system. The same is true for (2.2b) and 


(2.2d), but (2.2b) has additional constraints on the spin states and (2.2d) 


involves imaginary interaction coefficients. The form in (2.2c) is that of a 
vertex model with vertex weights given by cu,. 

We begin to present our results by first restricting to quantum systems 
with nearest neighbour interactions only. The extensions to longer-range 
interactions are detailed in Appendix [B] 


2.1 Nearest neighbour interactions 


Restricting (1.1) to nearest-neighbour interactions gives the well known one- 
dimensional quantum XY mode£] 


M 


h xv ,Y ( + k) • 


(2.3) 


3 = 1 


3 We can ignore boundary term effects since we are interested in the thermodynamic 
limit only. 
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where JJ = - (A jd+1 + ^B jd+l ), J y = - (A jd+1 - jB jd+1 ). This mapping 
is achieved by using Jordan-Wigner transformations: 


b) = ^ (m 2j+ 1 + im 2j ) b, = ^ (m 2j +i - im 2j ) 

= l (7 + *7) id (-O • = 5 (7 - *7) ff (-7), 

Z=1 Z=1 


(2.4) 


where 


m 2i+ i = a] JJ (-a]) , m 2j = a] JJ (-aj) , 
z=o 1=0 


or inversely as 


(2.5) 


3~ 1 


= 


im 2j m 2j+ 1 , o-J = m i2 +i JJ (-fm 2Z m 2i+ i), 


1=0 


3 - 1 


and crj = m 2i (-zm 2 i m 2 i+ i). 


z=o 


( 2 . 6 ) 


The rrijS are thus Hermitian and obey the anti-commutation relations {rrij, irtk} 
2 bjk- 


2.1.1 A class of classical Ising type models (2.2a) 


When we restrict to 7 = 1 and -£ 77+1 = Aj d+ 1 , (2.3) becomes a class of 
quantum Ising type models in a transverse magnetic field with site-dependent 
coupling parameters. Suzuki showed El that the partition function for such 
a system can be mappecQto that for a class of two-dimensional classical Ising 
models with Hamiltonian B rl given by 


n M 

n cl = - EE (Jj s j,p s j+ 1 ,p T J - 1 ) j ( 2 . 7 ) 

p = 1 3 =1 

with parameter relations 

(X,J V = \ logcoth 13 ' cl J? = ^ Jj, ( 2 . 8 ) 

2 n J n 

where /3 qu(cl) is the inverse temperature of the quantum (classical) system. 
4 Upto an overall constant. 
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Thus we have an equivalence between our class of quantum spin chains 
under these restrictions and a class of two-dimensional classical Ising models 
also with site-dependent coupling parameters in one direction and a constant 


coupling parameter in the other. From (2.8) we see that the magnetic field h 


driving the phase transition in the ground state of the quantum system plays 
the role of temperature /3 cl driving the finite temperature phase transition of 
the classical system. 

This mapping holds in the limit n —> oo, which would result in anisotropic 
couplings for the class of classical Ising models, unless we also take /3 qu —>■ oo. 
This therefore provides us with a connection between the ground state prop¬ 
erties of the class of quantum systems and the finite temperature properties 
of the classical. 

In this case we can also use this mapping to write the expectation value 
of any function / ({cr 2 }) with respect to the ground state of the class of 
quantum systems as 

(/(K})) qu = (/(M)L, (2.9) 

where (/ ({s})) cl is the finite temperature expectation of the corresponding 
function of classical spin variables with respect to the class of classical systems 


(2.7). 


Some examples of this are the spin correlation functions between two 
or more spins in the ground state of the class of quantum systems in the 
z direction, which can be interpreted as the equivalent correlator between 
classical spins in the same row of the corresponding class of classical systems 


(2.7); 


\ a j a j+r/ c 


— (■ Sj }P Sj + r tP ) 


cl ? 


n 

\j =i 


n 

\j =i 


'hV 


( 2 . 10 ) 


2.1.2 A class of classical Ising type models with additional con¬ 


straints on the spin states (2.2b) 


Similarly, the Trotter Suzuki mapping can be applied to the partition function 


for the XY model (2.3) in full generality. In this case we first order the terms 


in the partition function in the following way 


Z = lim Tr 




t> = X 


ftqu 


H 


ie 


Pojx 


x ftqu n / 


ie 


ftqi 


m 


( 2 , 11 ) 


j&a 


where H 1 ' = Jj<J^cr^ +l for fi G x,y, 77 J = | (cr| + cr| + i), and a denotes either 
a or b } which are the sets of odd and even integers respectively. 
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We then insert 2 n copies of the identity operator in the cr z basis; I s = 
|s p ) (s p | where |s p ) = |si p , s 2p ,..., sm p ) between each of the 2 n terms in 

pill ); 

Z = lim Tr I S1 V a I S2 V b ... I sa „_ 1 V 0 I a2n V6 


2 n 


lim 'y ' I I {sj p \ V a |sj iP _|_i) (syp+i| Vfe I) 

n .— 


s j,p P^ a 


The remaining matrix elements in (2.12) are given by 

M 


jea 


( 2 . 12 ) 


( s j,pl l s i,p+i) — ( s j,pi s i+i,pl M l s i,p+i> s j+i,p+i) j (2-13) 


where 


/ 


M = 




e cosh '^Bj j 


0 

0 

(4-4-1M | ft) 


V 


sinh 


^ 2/3qu7 


Bj 


0 

(niM|u> 


0 

(t4-|M|4-t> 


cosh [ ^AA sinh ( 


2/3, 

n 

atjMin) 

sinh ( ^^-Aj ) cosh 


n 3 

4t|M|4-t) 

2/3qu ^ 


(tt|M|U> 


sinh 


^2/3qu7 


Bj 


0 


0 

(U|M|U> 

cosh D, 




(2.14) 


It is then possible to write the terms (2.13) in exponential form as 


M 


( S j,pl l S i,P+l) — e 

jGa 


Pc\Hj,p 


(2.15) 


where 'Hj.p can be written as 


Bj, P — ^ (Jj s j,p s j,p+i + Jj s j,p s j+i,p + Jj s j+i,p s j,p+i 
+ H ( Sj )P + Sj+i, p ) + Cj), 


(2.16) 


or more symmetrically as 

Bj, P = — ( s j,p s j+i,p + s i,p+i s j+i,p+i) + Jj ( s j,p s j,p+ 1 + s j+i,p s j+i,p+i) 


+ Jj ( s j,p s j+i,p+i + s j,p+i s j+i,p) 

+ H ( Sj :P + Sj+i iP + Sj tP+ 1 + Sj+i,p+i) + Cj), 


(2.17) 
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where 


PciJ- 

/WJ 


P cl H 


log 


sinh ^ 7 Bj 


sinh ^ A- 

n J 


0c\Jj = log 


tanh 

n J 


tanh 


27A 


n J 


log coth '' qu Bj coth —^Aj 


n 


n 


31 


[X:\Cj — log sinh 


27& 


n 


n 


qu Aj sinh —AAL 


n 


Bj, 


(2.18) 


as long as we have the additional restriction that the four spins bordering 
each shaded square in Figure [l] obey 


Sj ,pSj +1 ,p Sj ,p+1 Sj +1 ,p+1 


1 . 


(2.19) 


This guarantees that each factor in the partition function is different from 
zero. 

Thus we obtain a partition function equivalent to that for a class of two- 
dimensional classical Ising type models on a M x 2 n lattice with classical 
Hamiltonian 7-L cl given by 


2 n M 


2 n M 


«•■= EE+EE K 


3,pi 


( 2 . 20 ) 


vpSa jda p£b j£b, 


(2.19) 


where T-Lj <p can have the form (2.16) or (2.17), with the additional constraint 


In this case we see that the classical spin variables at each site of the two- 
dimensional lattice only interact with other spins bordering the same shaded 
square, represented schematically in Figure [l] with an even number of these 
four interacting spins being spin up and down (from condition (2.19)). 
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Figure 1: Lattice representation of a class of classical systems equivalent to 
the general class of quantum systems (2.3). Spins only interact with other 
spins bordering the same shaded square. 


This mapping holds in the limit n —> oo, which would result in coupling 
parameters Jj 1 , Jf, H —> 0 and JJ —* oo unless we also take /3 qu —>■ oo. There¬ 
fore this again gives us a connection between the ground state properties of 
this class of quantum systems and the finite temperature properties of the 
classical systems. 


Again we have the same relationship between expectation values (2.9) 


and (2.10). 


2.2 A class of classical Ising type models with imaginary 


interaction coefficients (2.2d) 


Alternatively, lifting the restriction (2.19) we instead can obtain a class of 


classical systems described by a Hamiltonian containing imaginary interac¬ 
tion coefficients 


n M 

{Jj a j,p a j+i,p + Jj T i,p T j+i,p T j,p ( a j,p 

p= 1 i =1 


a iiP+1 )), (2.21) 
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with parameter relations given by 

Mi = —J*, A *JJ = —Jl A ,J = \ arctan-(2.22) 

J n J J n J 2 sink — h 

n 

To achieve this, we first apply the Trotter-Suzuki mapping to the quantum 
partition function divided in the following way 


Z = lim Tr 




(2.23) 


where this time = Ejli Jj a j a j+i f° r A 4 e x i V and ^ — Ej=i a j- 

Next insert n of each of the identity operators I ap = E<r, | a p ) (B p \ and 
I r = E r | t p ) (t p \, which are in the a x and cr y basis respectively, into (2.23) 
obtaining 

Z = lim Tr ... I T2r M '2 




lim TT<YII t p ) (t p | U 2 \& P +i) • 


(2.24) 


a j,p’ T j,p p— 4 


It is then possible to rewrite the remaining matrix elements in (2.24) as 
complex exponentials, 


{(Tr 


I Mi \f r ) (f r \ 4 |^) = e^(K"V«; + .)+«?) 

Fp> (-Tpl e 2n kp+l) 


, | uo>z 

X (crJ e 2™ n 


(2.25) 


_ (~<2 Mg ’ (2 (^tp+^tp+i)+^tp) + 2 ^S 3 =i T j,p( a j,p a j,p+ 1 ) 


where = £ i=1 J x a jtP a j+hp , H y p = E j= i J y Tj, P T j+ i, p , D = \ arctan ^ , 

n 

C — | cosh , and we have used 


(cr iiP |e a ^ |r iiP ) = ^cosh( 2 a)e*^ ctan (iiHEW) CT 2 >^>. 


(2.26) 


The classical system with Hamiltonain given by (2.21) can be depicted as 
in Figure [2j where the two types of classical spin variables <7j, n and r. 


can 


up 'j,p 

be visualised as each representing two-dimensional lattices on two separate 
planes, as shown on the left in Figure [2j The blue (thick solid) lines represent 
interactions with coefficients dictated by JJ, the red (thick dashed) lines 
by Jj, and the Jj coupling constants correspond to the green (thin solid) 
lines which connect these two lattice interaction planes. One can imagine 
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“unfolding” the three-dimensional interaction surface shown on the left in 
Figure [2] into the two-dimensional plane shown on the right. 



Lattice direction j -> 


n 1 

Y 

T 










| 

± 




1 

i 










.... 



1 i i 





fit 






-4 .Y- 




Trotter 


direction 


pi 


Figure 2: Lattice representation of a class of classical systems equivalent to 


the class of quantum systems (2.3). The picture on the left can be “unfolded” 
into that on the right. 


As in previous cases, this mapping holds in the limit n —> 00 , which 
would result in coupling parameters Jj 7 , Jj —» 00 and J —> unless we 
also take /3 qu —» 0. Therefore, it gives us a connection between the ground 
state properties of the class of quantum systems and the finite temperature 
properties of the classical ones. 

We can use this mapping to write the expectation value of any function 
/ ({cr 1 }) or / ({cr y }) with respect to the ground state of the class of quantum 
systems (|2.3) as 


(/(K})L = (/(W)) C 


</(K})L = (/(M)) c 


(2.27) 


where (/ ({cr})) cl and (/ ({t})) c1 are the finite temperature expectation values 
of the equivalent function of classical spin variables with respect to the class 
of classical systems (2.21 )j^] 

An example of this is the two-spin correlation function between spins in 


the ground state of the class of quantum systems (2.3) in the x and y direction 


which can be interpreted as the two-spin correlation function between spins 


5 Recall from the picture on the right in Figure [ 2 ] that the a and r represent alternate 
rows of the lattice. 
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in the same odd and even rows of the corresponding class of classical systems 
(2.7) respectively; 


( a j a j+r) ( <J i,p cr i+np) c i ’ ( 


°’i <T |+r-) qu ( T j,P T j+r,p) ci 


(2.28) 


2.3 A class of classical vertex models (2.2c) 


Another interpretation of the partition function obtained using the Trotter 
Suzuki mapping, following a similar method to that of |3], is that correspond¬ 
ing to a vertex model. 

This can be seen by applying the Trotter Suzuki mapping to the quantum 


partition function ordered as in (2.11) and inserting 2 n identity operators as 
in (2.12), with remaining matrix elements given once more by (2.13). This 


time, instead of writing them in exponential form as in (2.15), we interpret 


each matrix element as a weight corresponding to a different vertex configu¬ 
ration at every point (j,p) of the lattice; 


/3 M 

(sj )P \e n l s t,p+i) = J. ( s !>p, s j+i,p, s j,p+ 1, s j+i,p+i) • (2.29) 

j eo 


As such, the partition function can be thought of as corresponding to 
a class of two-dimensional classical vertex models on a (^ + nj x (4f + nj 
lattice as shown in Figure |4| with Mn vertices each with a weight 
cA (s jjP s j+ljP , Sj-p+i, s J+ i,p + i) given by one of the following 


w i(+l> +1> +1; +1) 
n 4 (-i,-i,-1,-1, 


i, +i, +i, —i) 

U3 J 5 ( + 1, — 1 , + 1 , — 1 ) 


— 1 , — 1 , + 1 , + 1 ) 


e n 


COsh ( —Fq I , 


n 


= e cosh ( ), 


= uj{(+1 , -1, -1, +1) = sinh 
= 1, +1, -1, +1) = cosh 


= Cc4(+1, +1, —1, —1) = sinh 



(2.30) 


thus leading to a class of 8-vertex models with the usual 8 possible respective 
vertex configurations as shown in Figure [3] 
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VV'VV'V'V'V'V' 

UJ\ UJ2 CJ3 CJ4 6J5 CJq (jJ'j (jJ% 

Figure 3: The 8 allowed vertex configurations. 

The values of these weights depend upon the column j = 1,... , M of 
the original lattice, thus each column has its own separate set of 8 weights, 
as represented by the different colours of the circles at the vertices in each 
column in Figure |4j 



Lattice direction j —>• 


Figure 4: Lattice representation demonstrating how configurations of spins 
on the dotted vertices (represented by arrows ti) give rise to arrow configu¬ 
rations about the solid vertices. 

Once again, this mapping holds in the limit n —» 00 , which would result in 
weights CU 3 , 6 ^ 4 , cuy, Wg —> 0 and weights uj\, uj 1 2 , Uq —» 1 unless we also take 
/3 qu —> 00 . It thus gives us a connection between the ground state properties 
of the class of quantum systems and the finite temperature properties of the 
corresponding classical systems. 


14 













2.4 Algebraic form for the classical partition function 

Finally one last form for the partition function can be obtained using the 


same method as in Section 2.1.2 such that the quantum partition function 


is mapped to one involving entries from matrices given by (2.14). This time 


however, instead of applying the extra constraint (2.19), we can write the 
partition function as 


n M 


n M 


z =isl e a nn+nn 


\p£a j£a p£b j&b/ 


2 A, 


[(1 s j,p s j+i,p) (1 d- s j,p s j,p+i ) cosh Ajj+i 


n 


A (1 s j,p s j+i,p) (1 s j,p s j,p+i) s i n h Ajj+i 


(2.31) 


n 


T (1 T Sj,pSj+i,p) (1 sinh 


2/3 qu 7 


qu ' B 


n 


+ (1 + Sj jP Sj iP + 1 ) (1 + i) e n hSj ’ p cosh 


3,3 + 1 

2/3qu7 


B 


n 


3 , 3 + 1J 


2.5 Longer range interactions 

The Trotter-Suzuki mapping can similarly be applied to the class of quan¬ 


tum systems (1.1) with longer range interactions, to obtain partition func¬ 


tions equivalent to classical systems with rather cumbersome descriptions, 
examples of which can be found in Appendix |Bj 


3 Method of coherent states 

To use the method of coherent states for spin operators S l = |cr®, we first 


apply the Jordan-Wigner transformations (2.4) once more to map the Hamil 


tonian (1.1) onto one involving Pauli operators a 1 , i G x,y,z ; 


n = o (( A + k + a j a \ 


1 <j<k<M 


( k-1 \ M ( 3 - 1 ) 

II ~ a i ) ~ h J2 aZ r 

1 = 3 +1 / 3 = 1 

We then construct a path integral expression for the quantum partition 


function for (3.1). First we divide the quantum partition function into n 
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pieces 


Z — Tr e 0Hqu = Tr 


g AT7^q U g A'r'Hqu g AT7iq U 


= V" 


(3,2) 


where Ar = - and V = e ArWqu . 

n 

Next we insert resolutions of the identity in the infinite set of spin coherent 
states |N) between each of the n factors in (3.2). The coherent states for 
spin operators, labeled by the continuous vector N in three-dimensions can 
be visualised as a classical spin (unit vector) pointing in direction N such 
that they have the property 


(N|S|N) = N. 


(3,3) 


They are constructed by applying a rotation operator to an initial state to 
obtain all the other states as described in [TS] such that we end up with 


(N| & |N) = -SN\ 


with N l s given by 


N = ( N x , N\N Z ), 0 < 6 < Tr, 0 < 0 < 2v r, 
= (sin 9 cos </>, sin 9 sin </>, cos 9). 


(3,4) 


(3.5) 


Inserting these states between the n factors in (3.2) and taking the limit 
77 , —y oo we obtain 


Z = 


»N(/ 8 ) 


' N(0) 


PN (r)e~^ dTH(N(T))-s B ^ 


(3,6) 


where % (N (r)) now has the form of a Hamiltonian corresponding to a two- 
dimensional classical system given by 

W(N(r)) = (N(r)|« qn |N(r)) 

= T N* (t) N% (t) 

1 <j<k<M 

k -1 M 

+ (Aj,k ~ lB jk ) N? (t) N% (t)) (-Af(r))-h^iV/(T), 

i U ' ( 3 - 7 ) 

= ^2 (^j,fccos(0j (r) - (r)) + Bjtf cos (<j>j (r) + (p k (r))) 

1 <j<k<M 

k -1 M 

x sin (0j (t)) sin (6> fe (r)) JJ (- cos (0; (r))) - h ^ cos (9j (r)). 

Z=j+1 j=l 
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The term 

S B = J^dT<N(T)|T| N(T)) (3 . 8) 

appears through the overlap between the coherent states at two infinitesi¬ 
mally separated steps At = r l+l — t*, and is purely imaginary. This is the 
appearance of the Berry phase in the action P3U3SI. Despite being imagi¬ 
nary, this term gives the correct equation of motion for spin systems [18] . 


4 Simultaneous diagonalisation of the quantum 
Hamiltonian and the transfer matrix 


This section presents a particular type of equivalence between one-dimensional 
quantum and two-dimensional classical models, established by commuting 
the quantum Hamiltonian with the transfer matrix of the classical system un¬ 
der certain parameter relations between the corresponding systems. Suzuki 
IEI used this method to prove an equivalence between the one-dimensional 
generalised quantum XY model and the two-dimensional Ising and dimer 
models under specific parameter restrictions between the two systems. In 
particular he proved that this equivalence holds when the quantum system 
is restricted to nearest neighbour or nearest and next nearest neighbour in¬ 
teractions. 

Here we extend the work of Suzuki ca, establishing this type of equiv¬ 
alence between the class of quantum spin chains (1.1) for all interaction 
lengths when restricted to U(N ) symmetry^ and the two-dimensional Ising 
and dimer models under certain restrictions amongst coupling parameters. 
For the Ising model we use both transfer matrices forming two separate sets 
of parameter relations under which the systems are equivalent. Where pos¬ 
sible we connect critical properties of the corresponding systems, providing 
a pathway with which to show that the critical properties of these classical 
systems are also influenced by symmetry. 

All discussions regarding the general class of quantum systems (1.1) in 
this section refer to the family corresponding to U(N ) symmetry only, in 
which case, we find that 

[ft q »,V d ] = 0, (4.1) 


under appropriate relationships amongst parameters of the quantum and 
classical systems, when V cl is the transfer matrix for either the two-dimensional 

6 which we see from Table 111 means the matrices Aj t k and Bj^ have Toeplitz structure. 
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Ising model with Hamiltonian given by 

N M 

"H ^ ^ ^ ^ (r/1 'S./J S? ) 1 .J "b 'I'l .] '^7.] \ 1 ) J (4-2) 

* j 

or the dimer model. 

A dimer is a rigid rod covering exactly two neighbouring vertices either 
vertically or horizontally. The model we refer to is one consisting of a square 
planar lattice with N rows and M columns, with an allowed configuration 
being when each of the NM vertices is covered exactly once such that 


2h + 2v = NM , 


(4.3) 


where h and v are the number of horizontal and vertical dimers respectively. 
The partition function is given by 


Z = V x h y 

allowed configs 



allowed configs 


(4.4) 


where x and y are the appropriate “activities” and a — A 

The transform used to diagonalise both of these classical systems as well 
as the class of quantum spin chains ( 1 . 1 ) can be written as 



_ i/K_ 

e 4 


Vm 


y, e 2 M q3 (bjU q + ibjVq'j , 
j 


Vq = 



(4.5) 


where the rj q s are the Fermi operators in which the systems are left in diagonal 
form. This diagonal form is given by Eq. (1.3) for the quantum system and 
for the transfer matrix for the Ising model 


V +(-) = (2sinh2AV 


T\ _y - ' 

2 g — 2^ q e q 


{vhq- 


(4.6) 


where K t — [3,J l and e q is the positive root o^] 

cosh e q = cosh 2/Tj* cosh 2 K 2 — sinh 2K* sinh 2 K 2 cos q. (4.7) 

7 The superscripts represent anticyclic and cyclic boundary conditions respectively. 

1 1 

8 This is for the symmetrisation V = V 1 2 V 2 Vf of the transfer matrix; the other 
possibility is with V' = VfVxVf where V x = (2sinh2 e~ K ^ <, V 2 = 
e K 2 <T '°i+i ) and tanh K* = e~ 2Ki . 
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The dimer model on a two-dimensional lattice was first solved by Kaste- 
leyn [7] via a combinatorial method reducing the problem to the evaluation 
of a Pfafhan. Lieb [IS] later formulated the dimer-monomer problem in terms 
of transfer matrices such that V cl = V 2 is left in the diagonal form given by 

v d = n (K (vivq + V-qV-q - + {1 + 2a? sin 2 q)^) , (4.8) 

0<q<n 


with 


A q = 2 a sin q (l + a 2 sin 2 


1 q )~ 2 • (4.9) 

For the class of quantum spin chains ( |l.f ) as well as each of these classical 
models we have that the ratio of terms in the transform (4.5) is given by 


2 V q Uq 



f Oq_ 

b q 

sin q 

cosh 2K J cos q— sinh 2K^ coth 2Ki 
< singTl—tanh2XJ : tanhi^2 cosq') 

„ tanh 2K? 

cos q— tanh K 2 tanh 2K J cos z q— sinh 2 ^2 
1 

a sin q 


for 1-L qu 
for V, 

for V', 
for V 2 d , 


(4.10) 


which we show in the following sections will provide us with relationships 
between parameters under which these classical systems are equivalent to 
the quantum systems. 


4.1 The Ising model with transfer matrix V 


We see from (4.10) that the Hamiltonian (1.1) commutes with the transfer 
matrix V if we require that 


a q sin q 

b q cosh 2 Kl cos q — sinh 2K\ coth 2 K 2 


(4.11) 


This provides us 


with the following relations between parameters under 
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which this equivalence liolcl^ 

sinh 2/1 f coth 2 /I 2 = — —^7 -- 

b(L) 


a(L) + 0(A) 


and — 7 ^- -A, = — coth 2/1 2 tanh K *, 


(4.13) 


a{L) + b(L) 


or inversely as 

cosh 2 K\ = , tanh 2 / 1 2 = - ^ \J\a{L)f - ( b(L )) 2 , (4.14) 


where 


«=o 


[V 1 ' 


a(L) = a(L)J2 L , HO = HO E 


/=0 


2 / + 1 


and 


z(o) = r. 


From (4.14) we see that this equivalence holds when 

a(L) 


/(L) 


> 1, a 


\L) < a 2 (L - 1) + 6 2 (L). 


(4.15) 

(4.16) 

(4.17) 


For L > 1, we also have the added restrictions on the parameters that 


r A —11 

2 


E & ( fc ) E 

fc=l Z=0 


k 

21 + 1 


L-l 


^(—l)*cos fc 2 * g + ^/(/) cos fc g = 0 , (4-18) 

fc=i 


2=1 


and 


L—2 


r + E a[k) cos k q + J2 a + E (97) E ( ■ ) ( _1 )* cosL = °> ^ 4 ‘ 19 ^ 

fc=l fe=l z=o ' ' i=l ' ' 

9 Here we have used the De Moivre’s Theorem and the binomial formula to rewrite the 
summations in a q and b q (1.5) as 


a„ = F - 


L If 1 


k —1 1=0 ' 2—0 


£ £ <-u 


'*"* cos fc_2i g, 


(4.12) 


bn = 


= tan q^2b(k) ^ 
fc=l z=o 


k 

21 + 1 


2 — 0 


£( ) (-l)" l cos fc --g. 
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which implies that all coefficients of cos* q for 0 < i < L in (4.18) and of 
cos* q for 0 < % < L — 1 in (4.19) are zerc 


10 


When only nearest neighbour interactions are present in (1.1) (L = 1), 
with a{k ) = b(k) = 0 for k ^ 1 we recover Suzuki’s result |15| . 

The critical properties of the class of quantum systems can be analysed 
from the dispersion relation (1.4), which under the above parameter restric¬ 
tions is given by 

|A 9 | = 2 P+1 |cos (i_1) q\ ((a(L) cos q + a(L — l)) 2 + b 2 (L) sin 2 g) 2 , (4.20) 

which is gapless for L > 1 for all parameter values. 

The critical temperature for the Ising model m is given by 


k r = k 2 , 


(4.21) 


which using (4.13) and (4.14) gives 


a(L) = ±a(L — 1). 


(4.22) 


This means that (4.20) becomes 
|A,| = 2 P+1 |a(L) cos (i_1) q\ ((cosq± l) 2 + 


sin 2 q 


(4.23) 


which is now gapless for all L > 1, and for L = 1, (4.22) is the well known 
critical value for the external field for the quantum XY model. 

The correlation function between two spins in the same row in the classical 
Ising model at finite temperature can also be written in terms of those in the 
ground state of the quantum model; 


= <*ol vr 


1 1 

2 rr x rr x V 2 
a j a j+r V 1 


^o) 


= ($o|VrV^ +r V 2 |4> 0 ) 


= ((vr^v? 


V 2 n x V 2 

v 1 a j+r v 1 


= cosh 2 K; (^a-J +r ) qu - sinh 2 K* <ajaj +r ) qu 
using the fact that (cjcrj+r) = ( a j a j+r) = Cf f° r r and 


T 0 = <f> 




(4.24) 


(4.25) 


10 For example setting the coefficient of (cos q) to zero implies that T = 

-ESS i] (-i) i «(2i)- 
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from (1.3), (4.1) and (4.6), where To is the eigenvector corresponding to the 


maximum eigenvalue of V and 4>o is the ground state eigenvector for the 


general class of quantum systems (1.1) (restricted to U(N ) symmetry). 


This implies that the correspondence between critical properties (i.e. cor¬ 
relation functions), is not limited to quantum systems with short range in¬ 
teractions (as Suzuki |15j found), but also holds for a more general class of 
quantum systems for a fixed relationship between the magnetic held and cou¬ 


pling parameters as dictated by (4.19) and (4.18), which we see from (4.20) 
results in a gapless system. 


4.2 The Ising model with transfer matrix V' 


From ( 4. 10| ) the Hamiltonian for the quantum spin chains (1.1) commutes 
with transfer matrix V' if we set 


CL n 


sin q (1 — tanh 2A'* tanh K 2 cos q) 


tanh 2 


(4.26) 


bq cos q — tanh K 2 tanh 2K\ cos 2 q s jnti 9R . 

This provides us with the following relations between parameters under 


which this equivalence holds when the class of quantum spin chains (1.1) has 
an interaction length L > 1; 


tanh 2K\ tanh K 2 = 
and a S L - X) = 1, 


b(L - 1) 


or inversely as 


sinlr K 2 = 


a(L) 


tanh2A'* = 


2(5*(L)) 

1 

~ a(L — 1) 


b(L) 

a(L) 

b(L-l) 

b(L-l) 

tanh 2K* X 

5-(L) 

sinli 2 K 2 

b(L-l) 

and 

i 



where 


a*(L) = a(L - 2) - a{L)J2 


1=0 


21 


l. 


From (4.28) we see that this equivalence holds when 
a(L) (2a*(L) + a(L )) < a 2 (L - 1). 


(4.27) 


(4.28) 


(4.29) 


(4.30) 
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When L > 2, we have further restrictions upon the parameters of the 


class of quantum systems ( 1 . 1 ), namely 


L—2 




Y K k ) cos fc q + Y b ( fc ) 


k= 1 
and 


k =1 


L-3 


2=0 


2 / + 1 


L—l 


i =1 


^ ( . ](-!)* cos fc “ 2 *g = 0 (4.31) 


r +Y kc ° sk q- Y a ^ ( 21 ) 1 cosk 2q 


k =1 
L 


k =1 


/=0 


(4.32) 


fc=i 


i=0 


i=2 


Y a ( fc ) I] ( 2 ll E ( J (- 1 ) 1 cosfc ' 24 9 = °- 


This implies that coefficients of cos* q for 0 < i < L — 1 in (4.31) and of cos* q 
for 0 < i < L — 2 in (4.32) are zero. 

Under these parameter restrictions, the dispersion relation is given by 

|A 9 1 = 2 P+1 |cos L-2 q\ ((cos q (a(L) cos q + a(L — 1)) + a*(L)) 2 

+ sin 2 q ( b(L ) cosg + b(L — 1)))U 

which is gapless for L > 2 for all parameter values. 


(4.33) 


The critical temperature for the Ising model (4.21) becomes 

— a(L - 1 ) =a*(L) + a(L), 


(4.34) 


using (4.27) and (4.28). 


Substituting (4.34) into (4.33) we obtain 

|Aq| = 2 P+1 |cos i_2 q\((a(L) cos q — a*(L )) 2 (cos q — l ) 2 

+ sin 2 q (b(L) cos q + b(L — 1))) 2 , 


(4.35) 


which we see is now gapless for all L >2 (for L = 2 this clearly corresponds 
to a critical value of T causing the energy gap to close). 

In this case we can once again write the correlation function for spins in 
the same row of the classical Ising model at finite temperature in terms of 
those in the ground state of the quantum model as 


/ X X \ _ l x X \ 

\ a j,k a j+r,k / Ia — \ a j a j+r/ > 


(4.36) 


where T ' 0 is the eigenvector corresponding to the maximum eigenvalue of V* 
and 

^0 = 4> 0 . (4.37) 
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Once more this implies that the correspondence between critical prop¬ 
erties such as correlation functions is not limited to quantum systems with 
short range interactions; it also holds for longer range interactions, for a 
fixed relationship between the magnetic held and coupling parameters which 
causes the systems to be gapless. 


4.3 The dimer model with transfer matrix Vy 


D 


In this case, when the class of quantum spin chains (1.1) has a maximum 


interaction length L > 1, it is possible to find relationships between parame¬ 
ters for which an equivalence is obtained between it and the two-dimensional 
dimer model. For details and examples see Appendix [Cj 

For a(k) = b(k ) = 0 for k > 2, we recover Suzuki’s result [T5] |. 
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Appendix A Symmetry classes 


Classical compact 
group 

Structure of matrices 

Aj,k {Bj,k) 

Matrix entries 

(MO,-.* 

U(N) 

ST 

1 

S' 

1 

e 

flj-fc, j,k>0 

0+(2N) 

a(j - k) + a(j + k ) 

g 0 if j = k = 0 

if 

either j = 0, k = l 
or j = l, k = 0 

9j—k T 9j+ki j, k > 0 

Sp(2N) 

a(j — k) — a(j + k + 2 ) 

9j—k jA — 0 

0 ± (2N + 1) 

a(j — h) T a(j + k + 1 ) 

9j—k ~F fiJ+fc+ij jA — 0 

0~{2N + 2) 

a(j — h) — a(j + k + 2 ) 

9j—k fiJ+fc+2> j,k ^ 0 


Table 1: The structure of functions a(j ) and b(j) dictating the entries of 
matrices A = A — 2hi and B = 7 B, which reflect the respective symmetry 
groups. The gis are the Fourier coefficients of the symbol g M ( 0 ) of M.m- 
Note that for all symmetry classes other than U(N), 7 = 0 and thus B = 0. 


Appendix B Longer range interactions 

B.0.1 Nearest and next nearest neighbour interactions 


The class of quantum systems (1.1) with nearest and next nearest neighbour 
interactions can be mappccf~ i ] onto 


M 


= _ £(yyy +1 + jjyy +1 - (yyy +2 + J'vp» +2 ) 


a 


3 +1 


3 =1 


(B.l) 


+ K)> 


where Jf = \ (A jtj+2 + 7 B jtj+2 ) and = \ (A j:j+2 - 7 B j>j+2 ) using the 
Jordan Wigner transformations (2.4). 


We apply the Trotter Suzuki mapping to the partition function for (B.l) 
with operators in the Hamiltonian ordered as 


Z = lim Tr 


Pq- 


U qjx P qu njz Pqu gjy Pqu gjU Pqu nj z 
rLa g 2n 1-1 g n 1-1 b g n n ' a g 2n rL g 


Pqu qjy Pqu 


Pqx. 




(B.2) 


11 Once again ignoring boundary term effects due to our interest in phenomena in the 
thermodynamic limit only. 
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where again a and b are the set of odd and even integers respectively, and 

^ = Efea (pj K^i + pip) - J'PIPP) and U z = h p =l 
for n e x,y and once more a denotes either a or b. 


For this model we need to insert 4n identity operators into (B.2). We use 


n in each of the a x and o v bases and 2n in the a z basis in the following way 


. Pqu n ~ lx _ ffqu n ~ /z ftin -ivtf _ P^injy pqn -i/z ^qu -i;: 

Z = limTr I^e n Ha I Sl e 2n H e n Hb I T1 e n K “I Sl e 2n H e n Hb 

n—>oc 


dm IIKW 






|s 2p _i) (s 2p _i| e'an Hb I Tp) (B.3) 


nl 


O j,p jTj,S j p—1 

Tp 


Bqu n ~,y 3 qu -1~IZ £qu -yi . , , 

r n \ e » Ha \s 2p ) (s 2p | e 2 ^ H e - Hb |cr p+ i)]. 


For this system it is then possible to rewrite the remaining matrix ele¬ 


ments in (B.3) in complex scalar exponential form by first writing 


/ —* | £_ TV ,Ji I —* \ / —* I ' 'J-J •* ' 'j_y t> | —* \ 

\<Tp\e n a |s 2p _i) (s 2 p_i| e 2 « n e » \t p ) 


n z 


-Hi 


W 2 p-i| e n Ha |s 2 p) (s 2 p| |o>) 

= e ^HS(p) e ^^(2p-l) e ^^(p) e ^^(p) e ^^n2p) e ^^(p) 
((Jp| S 2 p— l) (^2p—1| 7~p) (^pl ^2p) (^2p| ^p+l) 5 


^qu^z ftqu ^cr 


(B.4) 


where 7^1 (p) = J j ( a j,p a i+i,P + Vj+i^jp + J? + i<r j , P s j +i, P <? j+ 2,p), 

Pp) = Yp a i\ jV j C T j,P r j+i, P + r j+i,P T j+2, P )+ J, j J +i T j,P s j +i,P T j+2,p) and U z (p) = 
YljLi s j,p■ We can then evaluate the remaining matrix elements as 


(cTp| S 2 p—l) (^2p—1| Tp) (r 2 p— 1| S 2p ) ('S 2 p| ^p+l) 

1 M 

_ _ "I T p IjT ( ^j,2p—lH^j,2pH^j,p^j,2p — l ^J,p+1 ^2p~h^-j,p(^>j,2p Sj,2p— 1)) 

2 4M 11 
3=1 


(B.5) 


Thus we obtain a partition function with the same form as that cor¬ 
responding to a class of two-dimensional classical Ising type systems on a 
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M x An lattice with classical Hamiltonian T~i d given by 


— /3 c {H c i 


qu 

n 


p= i tea 


J? 


„ ( <J i,p Cr j+l,p + Cr jr'+l,p°'j+2,p) '^ + l Cr j,P S t + l,P 0 ’i+2,p 


+ 2 ( T t,P T i+l,P + r j+l,P r t+2,p) ^'+l r j,P S J+l,2p-l' r j+2,p'j 

jeb v y 

E (^" 2 ” ^ Tj,,pT - J+1 ’ p T t+ 1 ,P T j'+2,p) — ^+lH,P S j+l,2p r j+2,p^ 


,J: 


E ^ „ ( Cr j,p+l Cr j+l,p+l E °'j + l,p+l <J i+2,p+l) 


j£b 

~ , ^j+l CT t,P+ lS i+ 1 ,2p CT J+2,p+l)) 
n M 


EB 

p=i j=i 

M . 

ITT 


P<Ji m 


2 n A 


-T I Sy2p + 


/3 qu h , 27T 


E — s 


2 n 4 


’i.2p 


(B.6) 


E ^ ^ (BpE^p- 1 BiP+i , ^ 2 p E H,p (B;2p s i,2p-i))) E 4//A4 In 2. 


j=i 


A schematic representation of this model on a two-dimensional lattice is 
given in Figure [5j with a yellow border representing a unit cell which can 
be repeated in either direction. The horizontal and diagonal blue and red 
lines represent interaction coefficients J x , J' x and J y , J' y respectively and the 
imaginary interaction coefficients by the dotted green lines. There is also a 
complex magnetic field term |^/i± ^ applied to each site in every second 
row as represented by the black circles. 
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Lattice direction j —y 


Figure 5: Lattice representation of a class of classical systems equivalent to 


the class of quantum systems (1.1) restricted to nearest and next nearest 
neighbours. 


This mapping holds in the limit n —> oo, which would result in coupling 
parameters ^ J x , ^ J y , ^ J ' x , ^ J ,y ,^h -> 0 unless we also take /L u ->• 
oo. Therefore this gives us a connection between the ground state properties 
of the class of quantum systems and the finite temperature properties of the 
classical. 

Similarly to the nearest neighbour case, the partition function for this 
extended class of quantum systems can also be mapped to a class of classical 


vertex models (as we saw for the nearest neighbour case in Section 2.1) or 
a class of classical models with up to 6-spin interactions around a plaquette 
with some extra constraints applied to the model (as we saw for the nearest 


neighbour case in Section 2.1). We will not give the derivation of these as 


they are quite cumbersome and follow the same steps as outlined previously 
for the nearest neighbour cases, and instead include only the schematic rep¬ 
resentations of possible equivalent classical lattices. The interested reader 
can find the explicit computations in [5|. 

Firstly in Figure [6] we present a schematic representation of the latter of 
these two interpretations, a two-dimensional lattice of spins which interact 
with up to 6 other spins around the plaquettes shaded in grey. 
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Figure 6: Lattice representation of a class of classical systems equivalent to 
the class of quantum systems (1.1) restricted to nearest and next nearest 
neighbour interactions. The shaded areas indicate which particles interact 
together. 


To imagine what the corresponding vertex models would look like, picture 
a line protruding from the lattice points bordering the shaded region and 
meeting in the middle of it. A schematic representation of two possible 
options for this is shown in Figure [7} 


. 1 . „ . 


Figure 7: Possible vertex representations 
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B.l Long range interactions 

For completeness we include the description of a classical system obtained 
by applying the Trotter-Suzuki mapping to the partition function for the 


general class of quantum systems (1.1) without any restrictions. 


We can now apply the Trotter expansion (2.1) to the quantum partition 


function with operators in the Hamiltonian (3.1) ordered as 

M—l 


Z = lim Tr [ O 

n —>oo -1- 

3 = 1 


j qu o ')x Pqu 'jjx 

(e n ^ij+lg n n j,j+2 


Pqt- 


g n ^j,M g 2n(M— 1) 

M M-j 

P qu 


/?qu 


njz Pquqjy P<yi n,y ^qu-ijs /tx 

H e « H i’ M ...e " H j,j+ 2 e n (B.7) 

M-J-1 


= lim Tr [TT((TT e " ^J+ fe )e 271( ' M ^ 1)W ( "TT e ™ 

n —>-oo X X X X XX 

j=l /c=l /c =0 

where = J? k af < El mi’ (~^f) for ^hl/ and 77* = h YjLi °y 

For this model we need to insert 3 Mn identity operators, nM in each of 


/3qu 


the a x , a y and a z bases into (|B.7|) in the following way; 

M—l 


<M-j 


Z = lim Tr If] (I., ( II e ^ niM 1 


/8qu 


j =l V k =1 

'M-j -1 

x I [] e^I—lipr 

k =0 


.-1 M—l 

/ M—j \ 

[IIK^ii 

f T T P( i u a}a: \ 

| 1 e « Tw+* 

=0 i=l 

V fc=i / 


(B-8) 


ftqu 


,U Z 


'M-j -1 


e ” ~ j ' M k ) l T t+ip) ( T j+jp\ a j+jp+ 1 ))' 


b'+lp I 


k =0 


For this system it is then possible to rewrite the remaining matrix ele¬ 


ments in (B.8) in complex scalar exponential form by first writing 
tM-j 

<3qu 


(S i+ir \ II <= " ^ H ' !?**> 


. fc=l 
'M-j -1 


( s j+jpl ( e 71 fc J I T j+jp) ( T j+jp\ a j+jp+ 1 ) 

\ k =0 / 

= e ^Efci J («J, J+ fc(p)+«?j +fc (p)+Hpf=iy w ') (a j+j p| s i+ip ) 
( s j+jpl T j+jp) ( r j+ipl ^i+j'p+i) 


(B.9) 
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where 7^ ifc (p) = J2k= j+ i Jj^j^Kp II?=/+i (- s kp)> 

uikip) = e^+i J lk T o,P T k, P n?=/+i (- s ^, P ) and 

Up = h E ;=l cr| p . Finally evaluating the remaining terms as 




"Sp) (^p| Fp) (aI ffi>+1) 


2v/2 


M M 


JPJ g 4 ((1 ^j.pXl s j,p)+ r i,p(l s j,p) <T j,p+i T j,p) 
J=1 


(B-10) 


The partition function now has the same form as that of a class of two- 
dimensional classical Ising models on a Mx 3 Mn lattice with classical Hamil¬ 
tonian 'H c i given by 


n-1 M a M 

«. = EE<7 e ( Jj,k a j,j+jp (T k 

P =1 i =1 k=j +1 


j+jP 


k -1 


Jj,k T j,j+jp T k,j+jp) 


/=j+l 


A, 


'^ pJ+ 1 n(M-l) /Z 4 j Sjj+jp (B.ll) 


Z7T 


+ ^ (1 a j,j+jP + T j,j+JP + a jJ+jP S j,j+jP 
~ T jj+JP S jj+jP ~ °’j,j+jP+l T j,j+jp)) + n M dl 


2v^‘ 


A schematic representation of this class of classical systems on a two- 
dimensional lattice is given in Figure [8] where the blue and red lines represent 
interaction coefficients and J y - k respectively, the black lines are where 
they are both present and the imaginary interaction coefficients are given 
by the dotted green lines. The black circles also represent a complex field 
( n(M-i) h ~ t) ac tffig on each individual particle in every second row. 
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Figure 8: Lattice representation of a classical system equivalent to the general 
class of quantum systems. 

This mapping holds in the limit n —>• oo, which would result in coupling 
parameters^ J ^ k , ^-h —* 0 unless we also take /3 qu —>• oo. Therefore 

this gives us a connection between the ground state properties of the quantum 
system and the finite temperature properties of the classical. 


Appendix C Systems equivalent to the dimer 

model. 

We give here some explicit examples of relationships between parameters 


under which our general class of quantum spin chains (1.1) is equivalent to 
the two-dimensional classical dimer model using transfer matrix (4.8). 


When L — 1, from (4.10) we have 

1 # 


6(1) sin q 


a sm q T + a(l) cosg’ 
therefore it is not possible to establish an equivalence in this case 


(C.l) 
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• When L = 2 from (4.10) we have 


m 


-—, if T = —a(2), a(l) = 6(2) = 0, (C.2) 


a sing — 2a(2)sing’ 

thus the systems are equivalent under the parameter relations 

“ = W- r = -“( 2 >- “(i) = K2) = o. 


(C.3) 


When L — 3 from (4.10) we have 

1 6(1) — 6(3) + 6(2) cos g 


a sing 2 sing (a(2) + a(3) cos g) ’ 
if T = —a(2), a(l) = —a(3), 6(3) = 0, 

thus the systems are equivalent under the parameter relations 

2a(3) a(2) 6(1) - 6(3) 


(C.4) 


a = 


6(2) ’ a(3) 6(2) F ^ 

a(l) = —a(3), 6(3) =0. 


(C.5) 


Therefore we find that in general when L > 1, we can use (4.10) to prove 
that we have an equivalence if 


1 sin q 1 m Efi' 1 ( 2 ,h) Ho (!) (-Ip cos*- 2 *- 1 q 


081119 T + a(l)cosg + ^ =2 a(/c)^| ) (-l)'( 2 A ))sin 2 'gcos fe - 2i g 
We can write the sum in the denominator of (C.6[) as 


, (C.6) 


Y a ( 2 -?') + cos qY a ( 2 j + 


3 = 1 


3 = 1 

ml 

. 2 J j 


S in 2 g(5]o(2i) J]('')(-l)< s m 2 l<- 1 >, 
j= 1 *=1 




+ cosg Y a ^j + l )Y (-lysin 2 ^ 9 


l=i 


i=l 


fc=2 Z=1 


X)0WE(- 1 )'( 2 ,)™ :1< '' 1> 9C08‘- a ?). 


(C.7) 
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This gives us the following conditions 


r = -i>(2j) 


3 = 1 

a(l) = - a ( 2 J + 1) = 0. 

3 = 1 


(C.8) 


We can then rewrite the remaining terms in the denominator (C.7) as 


j =i 


sin 2 g(^a(2j) ^ cos“" q 


2—1 






i= 1 x / p =0 
3 / *-1 

3 ' i V+P r »r>a 2 P+ 1 , 


J=1 
f m — 1 1 


2=1 

3 


a ( 2 i + !) ( J " cos “" 1 '9 


p=0 


(C.9) 


a(2j + 1) ^ 


j=i 

[f] 


i=i 


2.7 + 1 

21 


i -i 


E(“ 1 ) P+ ' 


cos 


2(j-p-l)+l 


q) 


p =0 


' '2f ‘ ' 


E“®)E (a) Ei- 1 )”*' 


cos 


20-P-l) 


?)■ 


j=l 2=1 x ' p=0 

Finally we equate coefficients of matching powers of cos q in the numerator 
in (C.6) and denominator (C.9). For example, this demands that b(m) = 0. 
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